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Influence of quantum effects on the internal waves and the Rayleigh- Taylor insta- 
bility in plasma is investigated. It is shown that quantum pressure always stabilizes 
the RT instability. The problem is solved both in the limit of short-wavelength per- 
turbations and exactly for density profiles with layers of exponential stratification. 
In the case of stable stratification, quantum pressure modifies the dispersion rela- 
tion of the inertial waves. Because of the quantum effects, the internal waves may 
propagate in the transverse direction, which was impossible in the classical case. A 
specific form of pure quantum internal waves is obtained, which do not require any 
external gravitational field. 

Studies of quantum plasmas was initiated by Pines in the 1960's jl, 2], where the finite 
width of the electron wave function gives rise to dispersion, important in the high-density 
and/or low temperature regime. A number of quantum plasma studies has since appeared 
Ref. jsl, e.g.^ kinetic models of the quantum electrodynamical properties of nonthermal 
plasmas and covariant Wigner function descriptions of relativistic quantum plasmas 
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effects can give rise to completely new effects in plasma environments 
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unified picture. 

The above examples mainly focuses on oscillations in homogeneous quantum plasma 
backgrounds. However, quantum plasmas can in practice often involve nonuniform density 
profiles, which often develop in an real {e.g. in astrophysics) or effective {e.g. in inertial 
confined fusion) external gravitational field. In the classical case, stratified plasma in a 
gravitational field inevitably exhibits either inertial waves or the Rayleigh- Taylor (RT) in- 
stability depending on whether the stratification is stable or unstable |30|. The purpose of 
the present paper is to study influence of the quantum effects on the internal waves and the 
RT instability. Here we show that quantum pressure always stabilizes the RT instability. 
The stabilization has a meaning of an effective " quantum velocity" reducing the instability 
growth rate. In that sense the stabilization is similar to the RT stabilization by an abla- 



tion flow in inertial confined fusion 31 
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34| . We solve the problem in the limit of 



short-wavelength perturbations. We also find exact solutions to the RT stability problem 
for density profiles with layers of exponential stratification. In the case of stable stratifi- 
cation, quantum pressure modifies the dispersion relation of the inertial waves. Because of 
the quantum effects, the internal waves may propagate in the transverse direction, which 
was impossible in the classical case. Even more, we obtain specific form of pure quantum 
internal waves, which do not require any external gravitational field. The results could be 
of significance for astrophysical and intertial fusion plasmas. 

We start with equations for continuity and momentum transport in quantum plasmas in 



the MHD approximation taking into account gravitational field 35 
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where rrie, rrii are electron and ion masses. We consider small-scale effects taking the gravi- 
tational acceleration to be g = — (^z, where is a constant. The hydrodynamic equilibrium 
is determined by the balance of forces 
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However, the hydrodynamic equilibrium does not necessarily imply thermodynamic equilib- 
rium. In the incompressible limit of an essentially subsonic plasma dynamics it allows both 
stable and unstable density profiles. 

We consider small perturbations of the equilibrium according to (f{x,z,t) = (fo{z) + 
(^(z) exp{iujt + ikx), where ip denotes any of the fiuid variables, k is the perturbation wave 
number, and uj is the wave frequency. In the case of unstable stratification the frequency 
should be replaced by the instability growth rate a = iu. Here we are interested in the 
incompressible plasma dynamics with 

dp 



dt 
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for any Lagrangian plasma parcel. Incompressible fiow is typical both for the RT instability 
and the internal waves [s^; equation (jl]) holds for u/k <^ Cs, where Cg = [{dP/ dp)sY^'^ is 



the sound speed. Then, to first order in the perturbed quantities, the system (II])-(l2l) reads 
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We start with the limit of short wavelength perturbations, k/a ^ 1, where a = dhxp^/dz 
is the local inverse length scale of the density profile. In that case we can use the Wentzel- 
Kramers-Brillouin method j^] presenting perturbations in the form ip{x,z,t) = ipo{z) + 
ipexp{ikxX + ikzZ + iut). Then, in the short wavelength limit Eqs. (jH]) - (IE]) reduce to 
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Solving equations 
quantum plasma 



TTl) we find the dispersion relation for the internal waves in 
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where the classical frequency of internal waves is designated by u'^i = —{g/po)dpo/dz > 0. 
In the case of "horizontal" internal waves, = 0, k^j. = k, the dispersion relation takes the 
form similar to electrostatic and electromagnetic plasma waves 



(13) 



where plays the same role as the plasma frequency and 
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is the characteristic quantum speed. Even in the case of zero gravity with Ud = we obtain 
specific quantum internal waves with the dispersion relation 



U = Uqk^. 



(15) 



In the case of unstable stratification g • Vp < 0, or 

g dpo 



Po dz 



>0, 



(16) 



we obtain the RT instability instead of internal waves. The RT perturbations of short 
wavelength are strongly localized within the layer with the most steep density profile 



max a = dlnpo/dz, e.g. see 
to the mode with k. = 



31 



32], and the maximal instability growth rate corresponds 
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As we can see from (1171) . quantum effects always play a stabilizing role for the RT instability. 



In a sense, this stabilization is similar to the effects o 
in the context of the inertial confined fusion 31, 32 



ab 



33, 



ation flow and thermal conduction 



34| . Still, the ablation flow may 
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be destabilizing for perturbations of long wavelength /c/a -C 1 because of the additional 



32 



34j . On the contrary, quantum 



Darrieus-Landau instability of a deflagration front 
effects are always stabilizing. 

We can also find exact analytical solutions to Eqs. (j5]) - (IHl) for certain density pro- 
files. For example, exact solution is possible for plasma consisting of one or several layers 
with exponential stratification po oc exp(Q;2;), d{\n pq) / dz = a = const. After tedious but 
straightforward algebra, the system Eqs. ([5]) - ([8]) may be reduced to a single equation 

which has the same mathematical form as the respective classical result 39|] but with the 
"effective" growth rate modified by the quantum terms 

= ^' + uy. (19) 

Therefore, reproducing calculations for the RT instability growth rate in the classical case 
with we find the respective result for quantum plasma as 



cr 



' - - (20) 



Below we give some examples for the solution to Eq. fll8p . 

1) Bounded plasma layer. We consider density profile po oc exp{az) within the layer 
< z < L bounded by two rigid walls. The boundary conditions at the walls are Uz = 0. 
Then we have oc exp{pz) and Eq. (fT8|) reduces to 

jj^Uz + a/jUz ( — 1 J k'^Uz = (21) 



V^eflf 

with the solutions 
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/^i. = -f ± X + ^ -f^^ I (22) 
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The boundary condition -u^ = at z = 0, L is satisfied for 

/ 2 \ 1/2 

/^«,2 « ,2\ ™ 



k'-^-k^ = (23) 



4 / L 



which leads to the instability growth rate 

In the limit of short wavelength perturbations we recover our previous result Eq. fll7p : the 
mode number n has the physical meaning of the scaled wave number in the vertical direction, 
kz = nn/L. 
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2) A transitional layer. We consider a transitional layer of width L separating two 
uniform plasmas of different density pi and 9pi = piexp(aL) 

p = Pi for z < 0, 

p = Pi exp{az) for < z < L, 

p = Gpi for z> L. (25) 

In the uniform layers we obtain Uz oc exp{±kz) from Eq. (|T8l) . which leads to the boundary 
conditions 

^ = ±A;n, (26) 
dz 

at z = 0, L, respectively. Solving Eq. ffTSjl with boundary conditions Eq. fl26|) we come to a 
transcendental equation 

J^ = l + ^coth(/5L), (27) 



where 



/?=f^ + A;^-^pV''. (28) 



Solution to Eq. fl271) is unique and /? is real for 



A;< fc, = - ( Vl + 4/ln'e-2/ln0) . (29) 



For k > kc the value P becomes imaginary, and we obtain multiple solutions to Eq. (j2j 
similar to i^^, see Fig. [H Different branches of the plot at k > kc correspond to different 
mode numbers n = — Q, which show the number of zeros for the eigenfunction u^.. In the 
limit of short wavelength perturbations, /c/a ^ 1, A;L ^ 1, Eq. (!27j) goes over to (IMl) . In 
the opposite limit of long wavelength perturbations, fc/a-Cl, A;L-Cl, we find 

a' = l^.fc - Uy. (30) 

The first term in Eq. flSU]) is the RT instability growth rate at the discontinuity, while the 
second one stands for quantum stabilization. In principle, the term with quantum stabi- 
lization also contains small parameter k/a -C 1, and it should be omitted within the model 
of discontinuous density. However, the quantum term contains also another dimensionless 
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FIG. 1: Solution to Eq. ([27|) . o^fij ^fga^ versus the scaled wave number kja for the density drop 
= 8. Different branches of the plot at k > kc correspond to different mode numbers n = — 6. 

which plays the role of a quantum Froude number, and which may be either large or small 
depending on a particular problem. For this reason, the quantum term in Eq. (!30|) may 
be important even in the limit of long wavelength perturbations. Numerical solution to 
Eq. (1271) is shown in Fig. l2] for different values of the quantum Froude number. In the 
domain k > kc the solution is not unique; in Fig. l2] we presented the mode providing the 
maximal growth rate. As we can see from Fig. 12], in the case of the quantum Froude number 
above unity, Fvq > 1, stabihzation happens already for perturbations of short wavelengths, 
k/a < 1. 

The effect of quantum dispersion on the RT instability and internal waves in inhomoge- 
neous systems can be of relevance in astrophysical environments [24] and ultracold plasmas 



2l| . The stabilizing effect due to the collective version of Heisenberg's uncertaintly relation 
could even dominate the dynamics in very dense plasmas. However, for direct applications 
to such systems, more detailed calculations have to be made. Here we have indicated the 
principle dynamics of inhomogeneous fluids where such quantum effects can play a major 
role. Generalizations and more detailed calculations, incorporating e.g. magnetic pressure. 
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FIG. 2: Scaled instability growth rate oj ^fga^ Eq. ([TTt) . versus the scaled wave number kja for 
the density drop = 8 and the quantum Froude number Fr^ = 0—1. 

are left for future research. 
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